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+y> Abstract 

^ Previous results on multi-charged Baby Skyrmion solutions have pointed to a mod- 

_G ular structure, comprised of charge two rings and single charge one Skyrmions, which 

combine to form higher charged structures. In this paper we present numerical evi- 
dence which shows an alternative finite chain, multi-charged global energy minimum 
^ Baby Skymion solution. We then proceed from the infinite plane, to Baby Skyrmions 

on a cylinder and then a torus, to obtain the solutions of periodic Baby Skyrmions, of 
which periodic segments will correspond to sections of large charge Baby Skyrmions in 
the plane. 
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1 Introduction 



The Skyrme model is a (3 + 1) dimensional non-linear theory of pions and admits topo- 
logical soliton solutions which describe baryons. In this paper we are concerned with a lower 
dimension analogue, known as the Baby Skyrmion model pQ. This theory has the Lagrangian 
density 

C = l -d^4> ■ - ^(d^ x d v <f>) ■ (d^ x 0» - m 2 (l - 3 ), (1-1) 

where the field is a three component scalar, <p = (4>i,<f)2,4'3) and takes its value on the unit 
sphere 4> ■ <fi = 1. 



Theory (1.1) is a modified 0(3) cr-model which includes a fourth power derivative and a 



symmetry breaking potential. The fourth power derivative is commonly known as the Skyrme 



term. The Skyrme term coupled with the potential term gives rise to a scale in (1.1), and 
therefore by Derrick's theorem pj, the theory admits soliton solutions. The potential is 
analogous to the pion mass term in the (3 + 1) dimensional Skyrme model. The positive 
constant k pre-multiplying the Skyrme term, in combination with the mass, m, determines 
the size of the Baby Skyrmion as being proportional to ^Jn/m. For finite energy, <p has to 
be a constant at the spacial boundary, thus we require the condition 

lim cf) = (0,0,1). (1.2) 

IX^oo 

The physical space, is a one-point compactification IR 2 U {oo}, which is topologically equiva- 
lent to S 2 . This compactification of the physical space and the unitary nature of the target 
space, gives rise to the field configuration <fi at fixed time being a map 

: S 2 -> S 2 . (1.3) 

Such maps can be classified by the homotopy class ^(S* 2 ) = Z, which gives the theory a 
topological degree, otherwise known as a topological charge, B. The topological charge of 
the map 0, is found to be the pull back of an area form of the target S 2 , integrated over the 
physical space 

B = deg[0] = — / • x d 2 cf))d 2 x. (1.4) 
4tt J R 2 

The aim of this paper is to discuss the structure of static multi-charge solutions, focusing 
on the chain-like structure that we conjecture to be the minimum energy configuration of 
multi-charge Baby Skyrmions. Firstly, we shall outline the well known symmetry reduc- 
tion used to calculate charge one and two Baby Skyrmions, then in more detail discuss the 
numerical computation of higher charge structures, producing the minimum energy chain 
configuration. We then extend our investigation to include periodic spaces, firstly investi- 
gating infinite chains on R 1 x S 1 and then Baby Skymion crystal structures on the torus, 
T 2 . 
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2 Charge one and two 



Static Baby Skyrmions are energetic minima of the energy associated with (1.1), namely, 

.2 



E 



-d t ct> ■ d t ct> + y \d x <fr x d 2 cj>\ 2 + m\l- fa) ) <Rr. 



(2.1) 



The static field equation found from the variation of (2.1 ) is found to be a highly non-linear 
partial differential equation, 







didi4> + m 2 e 3 + \(f> 



(2.2) 



+ k 2 (didi4>(dj4>- dj4>) + di4>(didj4>- dj4>) — didj4>(di4>- dj4>) — dj<p(didi<p- dj4>)) 



where e3 = (0, 0, 1) and A is a suitable Lagrange multiplier to impose the condition 4>-(p = 1. 

The only way to solve for <fi is to use a numerical technique. Throughout this paper we 
used gradient flow to solve these partial differential equations. A very effective simplification 



of (2.1), and hence its corresponding field equation (2.2), is to take advantage of the un- 



derlying symmetry of an 0(2) rotation in the physical space, which can then be neutralised 
by an iso-rotation of the target space. This technique is extensively covered in [I], where 
it is shown that the field, 4>, can be re-expressed with axial symmetry in terms of a profile 
function f(r), and a polar angle 6 as 



(j)(x) = (sin / cos(B9 - x), sin / sm(B9 - x), cos /), 



(2.3) 



with a global phase x £ [0, 2tt). The profile function f(r) takes the values f(0) = ir and 
/(oo) = 0. Using this simplification the static energy ( 2.1[ ) becomes, 



E = 2tt 



\f + (i + K 2 f)^ sin 2 / + m 2 (l - cos /) ) rdr. 



(2.4) 



Varying (2.4) gives 



f" r + 



k 2 B 2 sin 2 f\ V k 2 B 2 2r ,k 2 B 2 , \ B 2 , f , 
' — — sin 2 / + f'-^r- sin(2/) - — sin(2/) ■ 



2r 



2r 



■m 2 r sin / = 0. 



For comparative consistency with [1J, we chose k = 1 and m 2 



(2.5) 

0.1. A gradient flow 

algorithm is then implemented using (2.5) to minimise (2.4). The minimum energies for 
B = 1 and B = 2 are found to be E\ = 19.66 and Ei = 36.90 respectively. These energies 
are within 0.03% agreement with the results of pp. The profile functions for these solutions 
are then used to generate the charge density plots in figure [1} 

Only B — 1,2 minimum energy Baby Skyrmion solutions can be calculated using the 
above method, due to the symmetry imposed by (2.3). For B > 2 this approach creates 



ring Baby Skymion structures, which are not stable to perturbations which break the axial 
symmetry. Therefore, to calculate B > 2 minimum energy Baby Skyrmions, full field two- 
dimensional gradient flow is required. 
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Figure 1: Two-dimensional plots of the topological charge density of Baby Skyrmions with 
B = 1 on the left and B = 2 on the right. 



3 Larger charge Baby Skyrmions B > 2 

As stated above, for B > 2 it is found [1] that the axially symmetric solutions are not 
stable to axial breaking perturbations in full field simulations. To proceed we solve the 



static field equation (2.2), using gradient flow. For initial conditions it is conducive to use 
stereographic coordinates 

W= (pl+l f\ (3.1) 

1+03 

A sum of n charge one Baby Skyrmions in stereographic coordinates, equally separated 
on a circle with arbitrary phase is used as the initial condition. This configuration is then 
numerically minimised using gradient flow on a 200 x 200 lattice with spacing of Ax = 0.2. 
As a comparison to validate the two dimensional minimisation, the minimum energies of 
the B = 1 and B = 2 are calculated to be E x = 19.65, E 2 = 36.90. These B = 1,2 
values differ slightly from the two dimensional results in [1],|4] and [3J, the latter results 
seem to differ by 1% compared to the greater accuracy axial solutions. Our results are in 
better agreement with the previous dimensionally reduced solution with E\ = 19.66, and 
have a closer agreement with the dimensionally reduced results of the collective literature. 
We believe our slightly greater accuracy is due to the use of fourth-order accurate finite 
differences in comparison to the second-order differences commonly used. 



3.1 5 = 3,4 

Minimisation of B = 3,4 initial conditions reveal energies of E 3 = 55.58, E± = 73.61, which 
are close to the previous literatures energies plus a 1% perceived discrepancy in their com- 
putations. We also gain similar charge density distributions as in pQ , shown in figure [2] 

The configurations for B = 1, 2, 3, 4 which we obtain are essentially identical to those of 
[I]. The first major difference between our results and the results contained in pp arises at 
B = 5. 
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Figure 2: Two-dimensional plots of the topological charge densities for the B = 3 Baby 
Skyrmion on the top, and the B = 4 Baby Skyrmion on the bottom. 

3.2 B = 5 

A configuration of five randomly orientated B = 1 Baby Skyrmions, equally spaced on a 
circle is used as the initial condition. This initial configuration gradient flows to a chain 
configuration, of energy E 5 = 92.02, differing from the 5 = 3 + 2 solution proposed in [Tj. 
By setting up a system of Baby Skyrmions as in figure |3j which possesses an underlying 
symmetry about the y-axis, we are able to minimise to a B = 3 + 2 solution, which has an 
energy of E = 92.41. 

Therefore we believe our chain solution is of a lower energy. A third justification of our 
results is shown in figure |4j Here it is clearly shown that the initial condition flows through 
a slightly perturbed 3 + 2 solution, and minimises to the chain configuration. This shows 
that the chain solution must be of a lower energy than the accepted 3 + 2 modular solution, 
which is a very long lived saddle point in the energetic landscape. Also, [5] has a B = 5 
static chain configuration for a potential term of the form m 2 (l — 03) ' 9 , but there seems to 
be no comment on this new configuration. 
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Figure 3: Two-dimensional plots of topological charge density, on the left five B = 1 Baby 
Skyrmions, used as initial conditions to create the 5 = 3 + 2 solution on the right. 




Figure 4: Two-dimensional plots of topological charge density of B = 5 Baby Skyrmions, 
with the previously accepted configuration of B = 3 + 2 on the left which then flows to the 
new B = 5 chain on the right. 
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3.3 B = 6 



The next configuration which differs from the common results of the literature is the B = 
6 case. Using a circle of six B = 1 Baby Skyrmions as the initial configuration, led to 
a minimum energy configuration similar to the B = 6 configuration presented in [I], a 
configuration of three B = 2 Baby Skyrmions located at the vertices of an equilateral triangle 
as in figure [5] We found the minimum energy of this triangular solution to be E§ = 110.30. 
This configuration does not match with the incremental structure common to all the previous 
lower charge configurations. This deviation from the chain form led us to minimise the energy 
of a row of six out of phase Baby Skyrmions, this produced a B = 6 chain of static energy 
Eq = 110.22. Due to the difference in energy being less than 0.1%, we believe the triangular 
structure to be a local minimum energy configuration. 




Figure 5: Two-dimensional plots of the topological charge densities of the previously ac- 
cepted B = 2 + 2 + 2 Baby Skyrmion triangular configuration on the left, and the new B = 6 
chain solution on the right. 



3.4 B > 6 

For charges 6 < B < 10 we systematically found the minimum energy solution to be a Baby 
Skyrmion chain. We also energetically minimised some larger charges, namely B = 20, 21, 
and again found chain solutions. 
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B 


E 


E/AtiB 


1 


19.65 


1.564 


2 


36.90 


1.468 


3 


55.58 


1.474 


4 


73.61 


1.464 


5 


92.02 


1.464 


6 


110.22 


1.462 


7 


128.55 


1.461 


8 


146.81 


1.460 


9 


165.11 


1.460 


10 


183.38 


1.459 


20 


366.20 


1.457 


21 


384.48 


1.457 



Table 1: Minimum energies of charged Baby Skyrmions 

Some features worth noting about the form of the chains, are that for large charge Baby 
Skyrmions the main central body of the chains appears to be made of uniform overlapped 
charge one Baby Skyrmions. Also the charge density tends to be less uniform nearer the 
ends of the chain, where upon they appear similar to highly perturbed charge two Baby 
Skyrmions. These B = 2 Skyrmions seem to 'cap' the ends of the chain, as in figure [6] for 
the B = 20. For small charges, the energy per charge does not decrease smoothly, as shown 
in figure [7} 




Figure 6: Two-dimensional plot of the topological charge density of the B = 20 Baby 
Skyrmion chain. 

Instead of a smooth decrease in energy for small B, the energy per-unit charge shows two 
trends, one for odd charge and one for even; these two trends then converge for increased B. 
This asymptotic energy per-unit charge and the periodic nature of the charge density of the 
central region of a long chain, as in figure [6j motivated the calculation of the energy per-unit 
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charge of an infinitely charged Baby Skyrmion chain, as described in the next section. 



Cg 1.465 

LU 



1.455 



even charge 
odd charge 
periodic chain charge 



B 

Figure 7: Plot of energy per-unit charge in units of 4ir, against the charge of the Baby 
Skyrmion chain. 



4 Baby Skyrmions on a cylinder. 

This section is similar to the work in [7], with the aim to calculate the minimum Baby 
Skyrmion energy per charge of an infinitely charged chain. As presented in [TJ, it can be 
shown that two well separated Baby Skyrmions, with a relative phase Ax, have an interaction 
energy which can be calculated by a dipole approximation to be 



2 2 

E int = K (mR) cos(Ax), 

7T 



(4.2) 



where Kq is the order zero modified Bessel function, p is a numerically found asymptotic de- 



cay constant and R is the large separation distance. (4.2) shows that the relative phase, Ax, 
of two Baby Skyrmions dictates the type of interaction. For Ax = 0, the Baby Skyrmions 
repel, and for Ax = tt they attract. Hence to create a stable chain of interacting Baby 
Skyrmions, we have to start with an anti-periodic chain of B — 1 Baby Skyrmions. To 
proceed, the same numerical procedure as outlined above is implemented, but instead of 
imposing (ft — e 3 a t the boundaries, we imposed periodic boundary conditions in the y- 
direction. We then vary the length of the periodic y-direction to minimise the energy as 
shown in figure [8j Due to the smaller lattice size than the infinite plane case, Aa; = 0.08 is 
used for a greater energy resolution. 

This shows a minimum energy for a cell of periodic length L = 8.56, for k = l,m = yCTl. 
The minimum energy for a Baby Skyrmion on a cylinder is found to be E m i n = 1.4549, this 
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Figure 8: Plot of Baby Skyrmion chain energy per unit charge, in units of 4rr, for a varying 
periodic cell length. 



energy is the dashed line plotted in Figure [7] which matches to the asymptotic of the energy- 
charge plot. 



5 Baby Skyrmion on a tours 

With the aim to discover the minimum energy per charge of a Baby Skyrmion crystal, 
periodic boundary conditions need to be imposed. Hence the physical space is now a 2-torus 
and so, 4> : T 2 -> S 2 . 

It is already known that the minimum energy configuration is that of a hexagonal lattice 
[8], where each hexagon hole effectively represents a charge of B = 1/2. For numerical 
simplicity it is easier to work on a fundamental rectangular torus of size L x y/3L, as in 0. 
This lattice has the geometry to contain eight hexagons, hence a B = 4 initial configuration 
is needed. To proceed we used the above algorithm to energetically minimise two B = 2 
Baby Skyrmions, on a L x \/3L bi-periodic lattice. The length L is then varied to find how 
the minimum energy varies with cell size, giving figure [9] 

The result is that the natural unit cell size for a Baby Skyrmion lattice is a cell of size 
L = 10.4, for k — l,m — \/oX Figure [9] also has the energy per unit charge data for the 
cylinder case where the periodic length is doubled to allow for a B = 4 unit cell. This shows 
that as the fundamental rectangle grows in size, locally the Baby Skyrmion behave as if it 
were on a cylinder. So the Skymion configuration for increasing fundamental rectangle, first 
starts off as a compressed crystal, then becomes an hexagonal crystal, then a deformed chain, 
then a energetically minimum chain and finally a single B = 4 Baby Skymion on an infinite 
plane. One of the main points of interest of the torus case, is that the minimum energy 
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Figure 9: Plot of Baby Skyrmion chain energy per-unit charge, in units of 47T, for varying 
periodic cell length. 



of a lattice E[ att = 1.4555 and the minimum B = 4 chain E c h ain = 1.4549 are well within 
0.05%. This difference in energy is smaller than the numerical accuracy hence it cannot be 
concluded which, if either, is the minimum energy solution, only that the two configurations 
have very comparable energies. 



6 Conclusion 

We have presented numerical solutions for multi-charge Baby Skyrmions which differ to the 
configurations in the literature. This discrepancy is believed to be due to the extra computing 
power now available, which allowed us to find lower energy solutions for B > 5. This has led 
to the proposal that Baby Skymion minimum energy configurations are chains, and we are 
confident in this result because of the smooth lowering of energy per charge for increasing 
charge configurations. Also, the energy per-unit charge asymptotes to the periodic chain 
energy per-unit charge case, showing the increasing stability as the charge of a chain grows. 
The planar hexagonal lattice and the periodic chains minimum energies are too similar to be 
able to identify the preferred configuration. This energetic similarity between the hexagonal 
lattice and the chain energies is similar to the negligible difference between the commonly 
accepted B = 2 + 2 + 2 solution and the B = 6 chain solution. The energetic similarity 
between the 5 = 2 + 2 + 2 and the B = 6 chain is an exact analogue to the full Skyrme 
model, where it is found that the B = 12 solution can be either a linear chain of three 
B = 4 Skyrmions, or a configuration of three B = 4 Skyrmions, each placed at a vertex of 
an equilateral triangle [TD]. These two configurations, again have similar energies. 

An interesting extension to this work would be to extend the numerical algorithms to 
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higher orders, with the aim to be able to confidently identify the minimum energy periodic 
charge configuration. Also it would be interesting to see if similar chain solutions are also 
minimum energy solutions to the full Skyrme model. 
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